Improving quantum gate fidelities by using a qubit 
to measure microwave pulse distortions 
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We present a new method for determining pulse imperfections and improving the single-gate 
fidelity in a superconducting qubit. By applying consecutive positive and negative n pulses, we 
amplify the qubit evolution due to microwave pulse distortion, which causes the qubit state to 
rotate around an axis perpendicular to the intended rotation axis. Measuring these rotations as a 
function of pulse period allows us to reconstruct the shape of the microwave pulse arriving at the 
sample. Using the extracted response to predistort the input signal, we are able to improve the 
pulse shapes and to reach an average single-qubit gate fidelity higher than 99.8%. 



A basic requirement for building a quantum informa- 
tion processor is the ability to perform fast and precise 
single- and two-qubit gate operations pQ. For qubits de- 
fined in superconducting circuits, much work has been 
done to improve the quality of both single-qubit [2H1] 
and two-qubit gate operations [5TU1|. Still, gate fideli- 
ties need to improve further to reach error rates small 
enough for practically implementing fault-tolerant quan- 
tum computing with error-correcting protocols [121 1 1 3] . 
In most qubit architectures, many single-qubit operations 
are implemented by applying short microwave pulses res- 
onant with the qubit transition frequency. The phase of 
the microwave pulse controls the rotation axis in the x-y 
plane of the Bloch sphere, whereas the pulse amplitude 
and duration set the rotation angle. A difficulty with 
this approach is that the single-qubit gate fidelity be- 
comes highly susceptible to any impedance mismatch in 
the microwave line between the signal generator and the 
qubit, since such imperfections lead to pulse distortions. 

Consider the microwave pulse shown in Fig.nTa), which 
initially has a Gaussian-shaped envelope A\\t) with a 
well-defined phase. When passing from the generator to 
the device, the pulse is distorted and acquires a quadra- 
ture component ^4q(£). The pulse was intended to per- 
form a rotation around the i-axis of the Bloch sphere [see 
Fig.jljb)], but the quadrature components present in the 
distorted pulse shape will change the rotation axis and 
generate an error in the final qubit state. The systematic 
errors due to the non-zero ^4q(£) are particularly prob- 
lematic for qubit control, since they will bring the qubit 
state out of the y-z plane expected from a pure rotation 
around the ir-axis. 

The distortion is described by the transfer function %, 
which is the frequency-domain representation of the sys- 
tem's impulse response h(t). If the transfer function is 
known, it is possible to correct pulse imperfections using 
digital signal processing techniques. By numerically ap- 




FIG. 1. (a) Distortion of a Gaussian-shaped microwave pulse. 
When passing through the transfer function rl, the pulse 
shape gets distorted and acquires a quadrature component 
Aq. (b) Bloch sphere describing the qubit dynamics in the 
rotation frame. The drive fields Ai and Aq generate rotations 
around the x- and y-axes, respectively. The distorted pulse 
shape in (a) will cause rotations around an axis non-parallel 
to x, thereby generating errors in the final qubit state. 



plying the inverse T-L~ l to the input signal x, the pulse is 
predistorted in precisely the right way to give the correct 
signal % [h^ 1 [x]\ — x at the device. The difficulty lies 
in obtaining %. Since superconducting qubits operate at 
millikclvin temperatures inside a dilution refrigerator, it 
is generally not possible to probe the signal arriving at 
the qubit directly with conventional instruments such as 
a network analyzer or a sampling oscilloscope. 

In this work, we take a different approach and use the 
qubit 's response to various pulses as a probe for determin- 
ing H JT3]. We have designed and implemented a pulse 
sequence aimed at obtaining the unwanted quadrature 
component of the signal arriving at the qubit. The se- 
quence consists of pairs of positive and negative n pulses 
around the i-axis; the reversing of the pulse direction 
amplifies the quadrature component of the signal and 
causes the qubit to slowly oscillate around the y-&xis. 
By measuring the rotation frequency for different pulse 
periods, we are able to extract the time dependence of 
those quadrature components. From the obtained signal 
we construct the inverse transfer function H , and use 
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it to numerically predistort the input signal. The result- 
ing pulse shapes give a significant reduction in the gate 
error rate, as determined in a randomized benchmarking 
experiment |15) . With optimized pulse shapes, we ex- 
tract an average gate fidelity higher than 99.8%, which, 
to our knowledge, is the highest gate fidelity reported so 
far for a superconducting qubit. 

We use a flux qubit [16] , consisting of a superconduct- 
ing loop interrupted by four Josephson junctions. Bi- 
ased at the optimal operation point, the qubit 's energy 
relaxation time is T\ — 12 /is, and the dephasing time 
is T 2 * = 2.5 jits (see Ref. [17] for a detailed device descrip- 
tion). The device is embedded in a SQUID, which is used 
as a sensitive magnetometer for qubit read-out [TH]. We 
implement the read-out by applying a short current pulse 
to the SQUID to determine its switching probability P S w 
When statically biasing the qubit loop at half a flux quan- 
tum < I'o/2 ($o — h/2e), the Hamiltonian becomes H = 
— | (ojqb <Tz + A{t) <7 X ), where = 5.4GHz is the 

qubit frequency and A{t) = A\{t) cos(wi) + A^{t) sin(wi) 
is the drive field. The drive is generated by applying an 
oscillatory flux $(t) to the qubit loop using an on-chip 
antenna, giving A(t) = 2Ip<&(t)/U, with I P = 180 nA be- 
ing the loop's persistent current. When driving the qubit 
resonantly (u> — to q b) and going to the rotating frame, we 
get 

H = ^(A l (t)a^-+A Q (t)af-) : (1) 

which is the Hamiltonian depicted in the Bloch sphere in 
Fig-gb). 

The microwave pulses are created by generating in- 
phase [Ar(£)] and quadrature L4q(£)] pulse envelopes 
using a Tektronix 5014 arbitrary waveform generator 
(AWG), and sending them to the internal IQ mixer 
of an Agilent 8267D microwave generator. We write 
the total transfer function from generator to qubit as 
T~L = Hext Hint, where H cx t refers to imperfections in the 
electronics and coaxial cables outside the cryostat, and 
■Hint describe signal distortion occurring inside the cryo- 
stat, for example from bonding wires or impedance mis- 
matches on the chip. To ensure that the pulses we send 
to the cryostat are initially free from distortion, we deter- 
mine Hext with a high-speed oscilloscope, and use 
to correct for imperfections in the AWG and in the IQ 
mixers [HI [20]. The setup allows us to create well-defined 
Gaussian-shaped microwave pulses with pulse widths as 
short as t pw = 2.5 ns [3T]. We define the Gaussian as 
A(t) — Aexp[— irt 2 /tp W ], so that the integrated area un- 
der the pulse (corresponding to the total rotation angle) 
equals At pw . The pulses are truncated at a total dura- 
tion of 3i pw . 

To extract information about Hint, we drive the qubit 
with consecutive pairs of positive and negative tt pulses 
in Ai(t), separated by the pulse period T. The sequence 
is depicted in Fig.^a), together with Bloch spheres de- 
scribing the qubit states at various points of the pulse 
sequence. Note that in Fig.[2|a), we show an example of 
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FIG. 2. (a) Pulse sequence used to probe quadrature com- 
ponents in the microwave pulses. By applying consecutive 
pairs of pulses with positive/negative amplitude to Ax, the 
qubit rotations due to the in-phase signal cancel out, leaving 
only contributions from the quadrature component Aq. (b) 
Bloch spheres depicting the evolution of the qubit during the 
pulse sequence in (a). The angle is the quadrature rotation 
acquired per tt pulse, (c) Qubit population after the pulse 
sequence in (a), measured vs pulse period T and total num- 
ber of pulses N, and projected onto the three axis x, y and 
z. The qubit undergoes slow rotations in the x-z plane due 
to the quadrature components in the microwave pulses. 



the drive pulses that reach the qubit, including a small 
AQ-distortion after each pulse to better illustrate how 
the sequence works. The signal we create at the generator 
does not have any quadrature components. Starting with 
the qubit in the ground state, we apply a 7r-pulse around 
x to take the qubit to |1) [step I in Fig.[2ja-b)]. Next, 
the ^Q-part, due to the pulse distortion, induces a small 
rotation 9 around y, bringing the qubit state slightly off 
the south pole (II). The negative tt pulse then takes the 
qubit back towards the north pole (III), but since this 
pulse is inverted, the following AQ-part rotates the state 
even further away from |0) (IV). After the first two pulses, 
the qubit has acquired a rotation of 29 around the y-axis 
(V). The sequence is then repeated, and for each pair of 
subsequent tt pulses the qubit rotates another 29. 

Figure [2^c) shows the qubit state after the pulse se- 
quence, measured versus the number of pulses and the 
pulse period T, and projected onto the three axes x, y, 
and z using additional tt/2 pulses to do state tomography 
before reading out the qubit's polarization [55]. There are 
clear oscillations in the x- and z-components, verifying 
that the qubit indeed rotates around the y-axis despite 
the pulses being applied to x. Note that the rotation fre- 
quency is relatively slow: it typically takes a few hundred 
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7r pulses to perform one full rotation around y. A striking 
feature of Fig.[2jc) is that the oscillation frequency varies 
with pulse period all the way up to T = 25 ns, much 
longer than the pulse width t pw — 2.5 ns. This indicates 
that the quadrature distortions persist for a substantial 
time after the pulse should have ended. 

To explain why the quadrature rotations depend on 
pulse period, we need to understand what happens when 
the 7r pulses start to overlap with the distortions of the 
previous 7r pulses. Let us start by assuming that the it 
pulses are instantaneous, and consider the qubit response 
to the static quadrature distortion shown in Fig.[3][a), 
where Aq(t) remains constant at Aq/2it = 0.4 MHz 
for 30ns after the 7r-pulse in A\(t) at t = 0. Figure 
[^b) shows the qubit quadrature rotation during the dis- 
tortion, plotted for different values of the pulse period 
T. If T is 30 ns or longer [black circles in Fig.[3]jb)], 
the qubit will continuously rotate in one direction dur- 
ing Aq(4), acquiring a total rotation per pulse of 9 = 

ft=o° nS A.Q(t)dt s» 4.3 deg. However, if the pulse period 
is only T — 15 ns [green squares in Fig.[3]jb)], the second 
7r pulse in A\ at t — 15 ns will reverse the direction of 
the ^Q-induced rotations of the first pulse, in the same 
way that a 7r pulse in a spin-echo experiment reverses the 
spin evolution due to low-frequency field fluctuations in 
its environment [23] • The rotation per pulse 9 acquired 
with pulse period T = 15 ns ends up being zero, since the 
rotations during the second half of An(i) exactly cancel 
out the rotations during the first half. For T = 10 ns [blue 
diamonds in Fig.[3jb)], there are two extra n pulses in A\ 
occurring during the distortions of the first pulse, and 
we end up with 9 — 1.4 deg. Note that we only consider 
the rotation due to the distortion of the first tt pulse; the 
total qubit rotation will be a sum of the rotations from 
all pulses. 

Having understood why 9 depends on pulse period 
T for a given Aq(4), we now ask if we can invert the 
problem: given a measurement of 9 as a function of T 
such as the black trace in Fig.[3jc), can we extract the 
pulse distortions v4q(£)? To simplify the problem, we 
discretize time in the smallest steps available with our 
AWG, At = 1/(1.2 GS/s) w 0.83 ns, and write A Q {t) as 

a vector Q = [Qx,Q2,- •■ ,Qn], with Q n — Aq(nAt). 
The rotations 9(T) in Fig.ga) are measured with the 
same time resolution, and in a similar fashion we write 
9(T) as 9 = [6>i, 6» 2 , • • • , 9 N ], 9 m = 9(mAt). Both vectors 
contain N = 30ns/At = 36 elements. We still assume 
the 7r pulses in A\ to be instantaneous, occurring with a 
period of m = T/ At in the discretized time. 

As explained previously, the ir pulses will act to peri- 
odically reverse the direction of the Q-induced rotations, 
and the total rotation angle 9 m generated by Q becomes 
a sum of forward and backward rotations, depending on 
the period of the n pulses: 
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FIG. 3. (a) Constant quadrature distortion Aq{£) used to 
illustrate how the rotation angle depends on pulse period, 
(b) Qubit quadrature rotation for the quadrature distortion 
shown in (a), calculated with the pulse sequence from Fig.[2|a) 
and plotted for a few different pulse periods T. (c) Quadra- 
ture rotation 6 acquired per 7r pulse. The black points are 
extracted from data similar to the measurement shown in 
Fig.J^c). The magenta points are the results using a pre- 
distorted pulse shape, aimed at minimizing the quadrature 
distortion, (d) Quadrature component ^4q appearing at the 
sample when applying a 2.5 ns wide Gaussian pulse A\ at 
the input of the experimental setup. The signal Aq is ex- 
tracted from the quadrature rotations shown in (c) (black 
data points). 



We can write Eq. (|8|) as a system of linear equations 9 = 
At M Q, where M is a matrix with elements being either 
1 or —1 depending on the direction of rotation [23]. By 
inverting the matrix, we get the quadrature distortions 
directly from the measured rotations 9: 

Q = M- 1 9/At. (3) 

In the experiment, the tt pulses have a finite width 
t pw = 2.5 ns. During the pulses, the qubit is strongly 
driven around the cc-axis, and the quadrature rotations 
due to Aq are suppressed to first order in Aq/Ai dur- 
ing those t pw /At = 3 discrete time steps when Aj >• 

Aq. The vectors therefore have to be limited to Q = 
[Qi,Q5, ■ ■ ■ ,Qn] and 9 = [#4,05, • • • ,9n], and the ma- 
trix M needs to be modified to include zeros at the po- 
sitions of the 7T pulses in A\ [25] . 

In Fig.[3]Jd), we show the extracted quadrature re- 
sponse Q, calculated using Eq. (J3J) and the rotation data 9 
from Fig.^c). For reference we also plot the shape and 
amplitude of the intended drive pulse A\, digitized at 
I n = Ai(nAf). The pulse has an amplitude of 200 MHz, 
giving a n rotation in t pw = 2.5 ns. The extracted 
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FIG. 4. (a) Qubit polarization as a function of the number of 
pulses in the RBM sequence, measured with and without pre- 
distortion. The pulses have width i pw = 3.5 ns, separated by 
the period T = 10.5 ns. We use the same number of sequences 
and randomizations as in Refs. [3] [15]. (b) Average error per 
pulse versus pulse width t pvj , measured with T = 3t pv . The 
predistorted pulses give a lower error rate for all pulse widths. 



quadrature response Q has considerably lower amplitude, 
but keeps oscillating for 25 ns after the main pulse ends. 

Next, we use the measured response shown in Fig.J^d) 
to determine the transfer function "Hi nt of the system 
|26j . With knowledge of H, we can calculate the inverse 
"Hr 1 and use it to predistort the input signal, with the 
aim of reducing the quadrature distortions. The magenta 
trace in Fig.[3jc) shows the quadrature rotations 9 for the 
same sequence of positive/negative pulses, but this time 
measured with a predistorted input signal. Compared to 
the black trace, 9 has been significantly reduced for all 
values of the pulse period T, thus validating our method 
and confirming that the pulse shown in Fig.[3]jd) actually 
corresponds to the signal appearing at the sample. We 
attribute the rotations still present after predistortion to 
errors due to oversimplifications in the linear model in 
Eq. ^ used to extract Q. It may be possible to get a 
better estimate for ^4q by calculating the qubit response 
using the full dynamics of the Hamiltonian in Eq. ([lj, 
but it would involve solving a system of 36 non-linear 
equations, which computationally is much harder than 
inverting the matrix M in Eq. ([3|. 

Note that there may also be pulse distortions appear- 
ing in the in-phase component A\. However, the con- 
secutive positive and negative pulses in the sequence of 
Fig.[2ja) will cancel the effect of any errors in the rota- 
tions around x, which is also confirmed in the exper- 
iment (the y-component in Fig.[2]jc) shows no oscilla- 
tions). This cancellation allows us to exclusively target 
the ^-distortions. 

Having determined a way to reduce quadrature distor- 
tions and improve the microwave pulse shapes, we pro- 
ceed to characterize the qubit gate fidelity in our system. 
A convenient way of testing single-qubit gates is to im- 
plement the randomized benchmarking protocol (RBM) 
|15j . where a random sequence of ir and 7r/2 pulses around 
the x- and y-axes are applied to the qubit. If the pulses 
are imperfect, the qubit will start to dephase as the pulse 
errors accumulate. Figure [4ja) shows examples of de- 
cay traces, where the three traces correspond to data 



measured with either full predistortion {% — 'Hcxt'Hint), 
with predistortion only for the room-temperature elec- 
tronics {% = "H ox t), or with no predistortion at all. The 
pulses with full predistortion give a significantly slower 
decay than those without; when fitting to an exponen- 
tial decay we find a decay constant of N = 537 ± 22 
pulses, corresponding to an average error per pulse of 
1/N = (0.186 ± 0.008)%, corresponding to a fidelity of 
0.9984. 

In Fig.Qb) we plot the average error per gate ver- 
sus pulse width t pw , with the pulse period set to T = 
3 i pw . The predistorted pulses perform better for all pulse 
widths, showing that the pulse shapes have improved and 
again confirming the validity of our method for deter- 
mining the transfer function Hint- The general trend is 
that the gate error is reduced for shorter pulses; this 
decreases the total time itotai — NT of the sequences, 
thereby reducing the errors due to loss of qubit coher- 
ence. The relevant coherence time during the RBM se- 
quence is a combination of T\, T2, and the coherence 
time during driven evolution; for simplicity we plot the 
expected error rate if the pulse errors were limited by 
Ti = 12 ns [dashed line in Fig.gb)]. This limit is a fac- 
tor of two lower than our best results, indicating that 
the predistorted pulses still contain some pulse imperfec- 
tions. We speculate that parts of the remaining errors are 
due to in-phase pulse distortions, which are not targeted 
with the method presented here. A similar scheme may 
be developed to investigate the in-phase errors indepen- 
dently. Another complication is that, in our system, J\ is 
strongly reduced when driving the qubit continuously at 
Rabi frequencies above 100 MHz, probably due to local 
heating. This may contribute to pulse errors for short 
pulses (where the drive amplitude A cx l/i pw becomes 
large). At high drive amplitudes the Bloch-Siegert shift 
will also start to introduce deviations from the rotating- 
wave approximation in Eq. ([!]) . 

To summarize, we have demonstrated a new technique 
of using a qubit to determine and correct microwave pulse 
imperfections, allowing us to generate single-qubit rota- 
tions with an average gate fidelity better than 99.8%. 
Even though there have been reports of superconducting 
qubits in 3D cavities with coherence times approaching 
100 /xs [27l [28], we note that we obtain a higher gate fi- 
delity in our system because we are able to create shorter 
pulses. By encoding the pulse imperfections into a slow 
rotation when applying many pulses, we are able to de- 
tect distortions on a nanosecond timescale without the 
need of a fast detector. This makes our method very 
general, and it can be applied to any quantum computing 
architecture where single-qubit gates are implemented by 
applying microwave pulses at the qubit frequency. 
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SUPPLEMENT SI: 

CORRECTING PULSE DISTORTIONS DUE TO INSTRUMENT IMPERFECTIONS 

The microwave pulses that implement qubit rotations are generated by creating in-phase (I) and quadrature (Q) 
pulse envelopes using a Tektronix 5014 arbitrary waveform generator (AWG), and sending them to the internal IQ 
mixer of an Agilent 8267D microwave generator. The AWG has a sampling rate of 1.2 GS/s and an analog bandwidth 
of 300 MHz. To ensure that the pulses we send to the sample are initially free from distortions, we use a high-speed 
oscilloscope (Tektronix DP071254B, with maximal sampling rate 40 GS/s and an analog bandwidth of 12.5 GHz) to 
sample the microwave pulses and extract the pulse envelopes. This enables us to infer the transfer function of the 
instrument, from which we use signal processing techniques to correct for imperfections in the AWG and in the IQ 
mixers of the microwave generator. 

In general, a transfer function hit) takes an input signal x{t) and converts it into the output signal y(t). In the 
time domain, the transfer function has the form of a convolution y(t) = h(t) * x(t); by going to Fourier space we get 

Y(f)=U(f)X(f). (4) 

or 

U(f)=Y(f)/X(f). (5) 

We can thus determine H(f) by measuring the response Y(f) for a known input signal X{f). If x(t) is a delta 
function, X(f) = 1, giving H(f) = Y(f), and the transfer function is equal to the impulse response of the system. For 
technical reasons, we found it easier to measure the step response instead of the impulse response, but the impulse 
response can easily be determined by taking the time derivative of the measured step response. To find the combined 
transfer function of the AWG and the IQ mixers, we used the following procedure (similar to Ref. [20J ): 

1. Generate a signal at the qubit frequency (/ mw = 5.4 GHz) with the microwave generator and send it to the 
mixer input. 

2. Apply a 5 /zs long square pulse to the mixer I channel, apply zeros to the mixer Q channel. The relatively long 
pulse duration is necessary to make sure the transfer functions captures the slow transients in the AWG output 
voltage. 

3. Sample the signal g(t) coming out of the mixer at 40 GS/s using the high-speed oscilloscope. 

4. Determine the in-phase L4i(i)] and quadrature [Aq(i)] components of the sampled signal. We do this by 
multiplying the sampled signal with a cosine (for Aj) or sine (for Aq) and integrate over one period T = l// mw 
to get the pulse envelopes. We have 

Ai{t) = ~f g{t') cos (27r/ mw t / ) dt' (6) 

2 r t+T 

A ^ t )^~Tj t S(*')sin(27r/ mw i') dt'. (7) 

5. Resample the envelopes A\(t) and Aq(£) at the sample rate of the AWG (in our case 1.2 GS/s). 

6. Convert the step responses A\{t) and Aq(£) to impulse responses by taking the numerical derivative. 

7. Create the complex impulse response function y(t) = (d/dt) (A\(t) + iAQ(t)). 

8. Calculate the Fourier transform of y(t) to get the transfer function H(f) = Y(f). Note that since y(t) is 
complex, H(f) does not generally have a symmetric spectrum. 

9. Create the inverse = 1/H(f). 

Once the transfer function has been determined, we use the inverse TL^ 1 to predistort the input signal. The system 
will then generate the desired output signal y(t) = h(t) * h^ 1 ^) * x(t) = x(t). Practically, we implement the following 
procedure: 

1. Generate (in MATLAB) the in-phase and quadrature L4q(£)] envelopes of the desired pulse sequence, 

sampled at the sample rate of the AWG. 
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2. Construct the complex input signal x(t) = A\{t) + iAq(t). 

3. Take the discrete-time Fourier transform of x(t) to get X(f). 

4. Predistort the signal by multiplying 'H~ 1 (f) with X(f). 

5. Return to the time domain by taking the inverse Fourier transform of T-L~ 1 (f)X(f). 

6. The real part of the resulting signal is the predistorted in-phase component, the imaginary part is the quadrature 
component. 

7. Generate the two signals in the AWG and send them to the I and Q ports of the mixer. 

Figure[5]demonstrates the improvement of the pulse shape after implementing the predistortion algorithm. In Fig.[5ja), 
we plot the extracted pulse envelopes of an intended 3 ns wide Gaussian pulse, created without any predistortion, 
while Fig. [5jc) shows the same pulse generated using predistortion. Figures [gb) and (d) are magnifications of low- 
amplitude regions of Figs^a) and (c), respectively. For the predistorted version, the ringing after the pulse and the 
quadrature components visible in Fig.[5^b) are practically gone. Also note that the pulse without predistortion in 
Fig.[5ja) does not reach up to the intended unity amplitude, which the predistorted curve does. 



SUPPLEMENT S2: 
DEFINING THE MATRIX M 



From Eq. (2) in the main text, we have 

m 2m 

n,„ = a/ X! c -- Y • Y °-- - Y 



3 m 



m+1 



2m+l 



3m+l 



(8) 



To see how this can be written in matrix form, we consider a specific example where the vectors Q and each have 
N = 10 elements. For m = 1 (where the period between ir pulses is one element), the total rotation 0\ becomes 



6>! = At [Q 1 -Q 2 + Q 3 -Q4 + Q5-Qe + Q7-Q8 + Q9- Qw] 



(9) 



and for m = 2 we get 



At [Qi + Q 2 - Q 3 - Qi + Qb + Qe ~ Q-i - Qh + Q9 + Qw] 



(10) 



(a) I 




20 30 



20 30 



Time [ns] 



Time [ns] 



W I 



0.5 



Ideal 
Gaussian 



10 

Time [ns] 



20 



0.04 
0.03 
0.02 
0.01 


-0.01 
-0.02 



30 




■ Aj 
.A Q 

. Ideal 
Gaussian 



10 20 
Time [ns] 



30 



FIG. 5. (a) Measured in-phase (I) and quadrature (Q) components of a 3 ns long Gaussian microwave pulse. The black dashed 
line shows the ideal Gaussian, (b) Magnification of the low-amplitude region of (a), (c-d) Same as panels (a-b), but here the 
pulse was generated using the predistortion algorithm. The pulse distortion that is strongly visible in (b) is essentially gone in 
panel (d). 



Written in this form, the matrix nature of Eq. pi) becomes more apparent. Writing out the full matrix, we have 



9 = At 
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SUPPLEMENT S3: 

MODIFYING THE MATRIX M TO INCORPORATE FINITE PULSE WIDTHS 



Equation (2) in the main text and the matrix M in Eq. (11) above both assume that the tt pulses reversing the 
rotation directions are instantaneous. This is not the case in the experiment, where we use Gaussian pulses with 
i pw = 2.5 ns. During the pulses, the qubit is strongly driven around the x-axis, and we can neglect the y-rotations 
from Aq at the t pw /At — 3 discrete time steps when A\ 3> Aq. This means that we can not extract any information 

about Q during the pulse, and the vectors will be limited to Q = [Q4, Q5, • • • , Qn] and 9 = [#4, 65, ■ ■ ■ , On]- We need 
to rewrite Eq. (2) in the main text as 



At 



.71=4 



m+4 



2m+4 



3m+4 



(12) 



The changes are most easily visualized by writing Eq. (12) in matrix form (again for the case N = 10): 



= At 



1 



0-1 



+10 

+1+10 

+1 +1 +1 

+1 +1 +1 +1 

+1 +1 +1 +1 







-1 
-1 





+1 



Q 



(13) 



Here, we have removed the entries for n,m = 1,2,3 (since the vectors are three elements shorter). Also, we have 
include zeros at the positions where the tt pulses occur in A\ (since during the 7r pulses, A\ 3> Aq 1 and the quadrature 
component is very inefficient in driving rotations around y). 



SUPPLEMENT S4: 

PREDISTORTING USING THE RESPONSE EXTRACTED FROM THE QUBIT QUADRATURE 

ROTATIONS 

The procedure for predistorting signals based on the qubit response is very similar to the sequence described in 
section SI, with some slight differences in how we determine the transfer function H . After having calculated the 
quadrature component Q as described in the main text, we proceed as follows: 

1. Construct the output signal y(t) = I + iQ from the two vectors / and Q. Here, Q is the extracted quadrature 
component, while we set / to be an ideal Gaussian with pulse width t pw — 2.5 ns and amplitude Ajlix — 200 MHz 
(our method does not give any information about in-phase distortions). 

2. Construct Y(f) by taking the discrete-time Fourier transform of the signal y(t). 
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3. In this case, the input signal x(t) is not an ideal delta or step function, but a Gaussian envelope with pulse 
width t pw — 2.5 ns and amplitude A/2ir = 200 MHz. Create x(t) by sampling the Gaussian at the same sample 
rate as the input y(t). The complex part of x(t) is zero. 

4. Construct X(f) by calculting the discrete-time Fourier transform of the signal x(t). 

5. Calculate the transfer function H(f) = Y(f)/X(f), as given by Eq. 

6. Create the inverse %~ l {f) = 

Once we have determined the transfer function and its inverse, we follow the same sequence as described in section 
SI for creating the predistorted signal. 

I 



